;)x-broadening and production processes 
versus dipole/quadrupole amplitudes at 
next-to-leading order 



A. H. Mueller 

Department of Physics, Columbia University, New York, USA 

S. Munier 

Centre de physique theorique, Ecole Poly technique, CNRS, Palaiseau, France 



Abstract 

Through the systematic inspection of graphs in the framework of hghtcone perturba- 
tion theory, we demonstrate that an identity between the evolution of p_L-broadening 
ampUtudes with the energy and the evolution of forward scattering amplitudes of 
color dipoles off nuclei holds at next-to-leading order accuracy. In the general case, 
the relation is not a graph-by-graph correspondence, neither does it hold strictly 
speaking for definite values of the momenta: Instead, it relates classes of graphs of 
similar topologies, and in some cases, the matching requires an analytical continu- 
ation in the appropriate longitudinal momentum variable. We check that the same 
kind of relation is also true at next-to-leading order between amplitudes for the 
production of dijets and quadrupole forward amplitudes. 



1 Introduction 

Observables such as dijet correlations [T|2] or other production processes in 
proton-nucleus scattering are outstanding probes of the high-density regime 
of QCD, and are likely to become one of the major focusses of the QCD 
community in the next few years. 

Preliminary phenomenological works [3] have shown that saturation models 
have the ability to describe successfully the present data. In order to confirm 
these results and provide accurate predictions for the Large Hadron Collider, 
firm theoretical foundations for the calculation of these production processes 
are now needed. 
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One outstanding problem in this program is the energy (or equivalently, the 
Bjorken-x) dependence of the considered observables. The formahsm to com- 
pute the energy evolution of forward amplitudes has reached a certain degree 
of sophistication: The evolution is known to be given by the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation [Hl5|6] . whose kernel is established at next- 
to-leading order in Inl/x [7|l8] see e.g. [9|T0] for references and for the state 
of the art of the phenomenology of deep-inelastic scattering. But production 
processes have not received the same attention so far. 

Interesting relations between production processes and forward amplitudes 
were conjectured some time ago (see e.g. Ref. [llj) and were proved at leading 
order [12j in particular cases. In this paper, we establish the relation between 
Px-broadening cross sections and dipole forward amplitudes [13] at next-to- 
leading order. 

To this aim, we examine essentially all possible lightcone perturbation theory 
diagrams [H] which contribute to these two processes, in a similar manner 
as was done to better understand the leading-order evolution of dipoles |T5] . 
and compare them. More precisely, our work consists in the discussion of the 
corrections to the wave functions up to order a^lnl/x and in the large num- 
ber of colors (Nc) limit. It is then clear that the arguments may be promoted, 
by iteration, to a fully resummed next-to- leading order result, namely to all 
diagrams up to order as{as Inl/x)^ for arbitrary k (with the further approxi- 
mation that quarks are left out of the evolution): Indeed, the graphs we review 
actually form the kernel of the energy evolution equation, which is known to 
be the BFKL equation on the dipole side [13] . 

Many graphs contribute at next-to-leading order in lightcone perturbation 
theory (although the large- A^^c limit reduces significantly their number by elim- 
inating the nonplanar ones). This is a drawback of time-ordered perturbation 
theory with respect to a covariant formalism. However, this formalism is the 
simplest and the most natural one to formulate the BFKL evolution, and thus, 
we believe that it is the most adequate for our purpose. 

The outline goes as follows. Section [2] presents a general discussion of the 
relations between production cross sections and scattering amplitudes. We 
then specialize to p^-broadening cross sections versus dipole amplitudes: Sec. [3] 
is devoted to a comprehensive review of their relation when a leading-order 
quantum correction is added, while the main new results of the paper on the 
relation at next-to-leading order are presented in Sec. |4} We briefly explain 
how these results may go over to other processes in Sec. |5| and conclude in 
Sec. El 
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2 Production cross sections and scattering amplitudes 

In this section, we shall discuss in a general way how certain production cross 
sections may be related to particular scattering amplitudes. There is of course 
one well-known case where this happens. Namely, the total cross section for 
the collision of two particles, A and B, is given by the imaginary part of 
the forward elastic scattering amplitude for particles A and B. This is the 
optical theorem and is valid in any quantum field theory. The relations that 
we are investigating in this paper, between cross sections for jet p^-broadening 
and jet production at high energy with dipole and quadrupole amplitudes, 
are specific to gauge theories as we shall see a little later on. There are also 
well-known relationships between inclusive particle production and particular 
discontinuities of higher-particle amplitudes. For example the cross section for 
A + B C + anything is, in general, given by a particular discontinuity of 
the elastic forward scattering amplitude A + B + C A + B + C This type 
of relation is true for the processes we are interested in here, but it is not this 
type of relationship that we are after. We are looking for relations between 
cross sections and the values of scattering amplitudes and in general there is no 
simple relationship between particular discontinuities of scattering amplitudes 
and the values of the amplitudes. To see in more detail the relationships we 
are investigating, it is perhaps useful to work through a simple example and 
then turn to the more general issues. 



2.1 p± broadening in the McLerran-Venugopalan model 

The simplest example of the type of relationship we are investigating is well 
illustrated in a formula which gives the p^-broadening of a quark (jet) as it 
passes through nuclear matter in terms of the scattering amplitude for a quark- 
antiquark dipole with the same nuclear matter. In the McLerran-Venugopalan 
model [16] the nucleons in a nucleus are uncorrected and the gluon distribu- 
tion of a nucleon xG{x, Q^) is evaluated at lowest order in DGLAP evolution 
with no x-evolution at all. 

Suppose a high-energy quark passes through a length L of nuclear matter. 
If the quark enters the nuclear matter with no transverse momentum {pj_ = 
0), and exits the nuclear matter with transverse momentum p^, then the 
distribution of transverse momenta is given by 




(1) 
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(a) 





Fig. 1. The three different types of interactions of a quark and a nucleus: inelastic 
(a), elastic in the amplitude (b) and elastic in the complex-conjugate amplitude (c). 
This illustrates the 3 terms in the r.h.s. of Eq. ([5]). 



where 



with 



(2) 
(3) 

The state 1^4) is the ground state of the nucleus while In) represents any final 
state after the quark has passed through the nucleus. The quark is represented 
by the Wilson line V which we take to move along the x_ = lightcone. y4_ 
and V are color matrices with the trace in Eq. ^ referring to the matrix 
indices in V^^ and Vqj^. The time-ordering (path-ordering) in Eq. ([s]) is such 
that A's at later times are placed to the left of A's at earlier times. At the 
moment S{x±) is simply defined by (|2|, however, we shall soon see that S is 
in fact the S'-matrix for dipole scattering on the nucleus. 

Because the nucleons in the nucleus are uncorrelated in the McLerran-Venugopalan 
model one has, for a sufficiently large nucleus. 



S{x±) 



oSiix±_) 



(4) 



where Si{x±) is the term in Eq. ^ involving scattering with only one nucleon 
in the nucleus. (The scattering in ^ occurs at a definite impact parameter 
although, for simplicity, we suppress the dependence.) Expanding Eq. ^ to 
order g"^ and using completeness of the states n one gets 



'S'i(a;_L) = 9^ I dx+dx'^{A\—ii\A^{{^^,x'^)A^{x±,x, 



1^ 



f(A_(0^,x+)A_(0^,x'+)) |A) (5) 



The three terms in the brackets {} correspond to an inelastic collision with a 
nucleon in the nucleus, an elastic scattering in the amplitude and an elastic 
scattering in the complex conjugate amplitude respectively as illustrated in 
Fig. [Ij The state N in the fi gure is that of a single nucleon in the nucleus 
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while the x+ and x', integrals in Eq. (jsj) go over the extent of the nucleus. 



In the McLerran-Venugopalan model the quark-nucleon scattering amplitude 
is purely absorptive so that 

{A\^tTA_ (xi_L, Xi+)A_ {X2±,X2+)\A) 

= {A\^tTT[A_{xi^,x,+)A_{x2±,X2+)] \A) (6) 
matching onto conventional notations when x^ is small 

(A I ^tr I A_ (Ox , O A_ (xx , a;+) 

1 1 1 

- ^A_{x^,x'^)A_{x^,x+) - 2^-(0±,x'+)A_(0x,x+) MA) 

71 1 

= -^xIp^G{x,1/xI)-^5{x+-x'^) (7) 

S^ = (8) 

Ql = ^^^LpxG{x,l/xl) (9) 

and where Qs is the quark saturation momentum and p the density of nucleons 
in the nucleus. L is the path length of nuclear matter at the impact parameter 
at which the quark passes through the nucleus. 

Using Eq. ^ in Eq. Q gives the relation 

S{x±) = e-"'^«'/^ (10) 

which is recognized as the elastic S'-matrix for a quark-antiquark dipole passing 
through the nucleus. Using Eq. (10) in Eq. ([T| tells us that the distribution 



so that 



with 



of transverse momenta that a quark picks up in passing through a nucleus is 
simply related, by Fourier transform, to the elastic dipole-nucleus S'-matrix 
element. If one integrates dN/(Pp over all transverse momenta the result must 
be 1, 

since the quark must pick up some value of the transverse momentum. In 
Eq. ([T]) this relationship corresponds to S'(0_l) = 1. From Eq. ( p|) 5 '(0_l) = 1 
corresponds to the unitarity of the operator Vq^ while from Eq. (|10[) it corre- 
sponds to color transparency. But color transparency is a property of gauge 
theories. Thus the relationships that we are investigating will not be true in a 
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general field theory but are specific to gauge theories. Finally, we have obtained 
Eq. (10) assuming that the dipole-nucleus scattering amplitude is purely ab- 
sorptive. While this is natural in the McLerran-Venugopalan model, where the 
cross-section has no energy dependence, the assumption is not necessary as 
the second and third terms on the right-hand side of Eq. (|5| combine to give 
an absorptive part in contrast to the individual terms which are time-ordered. 



2.2 p±-broadening more generally 



We now turn to a more general discussion of the possible relationship be- 
tween px-broadening and dipole scattering. We shall restrict our discussion to 
the case where gluons giving quantum contributions to p^-broadening have a 
longitudinal momentum much less than that of the quark, thus allowing the 
quark motion to be given by a Wilson line with a straight line trajectory. This 
restriction should not limit our ability to investigate and compare the energy 
evolution of p^-broadening with that of dipole scattering. A complete study 
of px-broadening vs dipole scattering would require an understanding of how 
to include the high-momentum gluons into (evolution-independent) impact 
factors and is beyond our present aim. 

Equation ([T]) is general but, when quantum gluon corrections are included, 
Eq. ^ should be changed to 

(12) 

where Cj is the QCD interaction Lagrangian, say in the lightcone gauge. S{x±) 
now has an energy dependence which, for simplicity, we suppress. Using com- 
pleteness of the states |n), which now include arbitrary numbers of quarks and 



gluons, in addition to nuclear break-up, one can write Eq. (12) as 



S{x±) = ^tr(A|f (V^^e-'l'^'y^'^y^) T (v^^e' I '"'^^'^y^) \A). (13) 
On the other hand, the S'-matrix for dipole-nucleus scattering is 

S{x^) = i^tr(A|T «K^e^/^'^^^(^)) \A). (14) 



Equation ( 14 ) is given as a conventional time-ordered product while Eq. ( 13 ) is 



a time-ordered product in a Keldysh-Schwinger [T7fT8] sense where the vertices 



in the T(), the second factor in the matrix element in (13), occur in the upper 
{xi+) part of the time contour illustrated in Fig. 2l while the vertices in T(), 
the first factor in the matrix element in Eq. ([l3), occur on the lower (x2+) 
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Fig. 2. Keldysh-Schwinger contour. 



A 




Fig. 3. The Feynman amplitude described by Eq. (15). 



part of the time contour p~|Feynman hnes connecting points on the Xi+-contour 
have ordinary ie factors, hnes connecting points on the X2+ part of the contour 
have the sign of it changed, corresponding to the anti-time-ordering indicated 



in Eq. (13). Feynman hnes connecting points on the Xi+ part of the contour 



to points on the X2+ part of the contour are put on shelL 



S and S given by Eq. (13) and ( |14[ ) respectively will not in general be the 
same. However, Kovchegov and Tuchin have shown [12] (see Ref. [20] for a 
broader review, see also Ref. [H]) that the x-evolution (energy evolution) 
of production cross sections is the same as dipole scattering at the leading 
logarithmic level. Our goal in this paper is to show that this equality persists 
through next-to-leading order in x-evolution. 



Finally, we note that ( 13 ) can be obtained as a particular discontinuity of an 



analytically continued Feynman amplitude. To see this we must keep track 
also of longitudinal momenta of the initial and final quarks as illustrated for 
the amplitude A in Fig. |3| 

P is the momentum of the target nucleus with p and p' the momenta of the 
quark. One can view A as having the following dependence: 



A = A({p+ p)^ {p-p'+ p)\ {p^ - p'^f 



Then 



S 



2i 



-disc 



(15) 
(16) 



as illustrated in Fig. |4} 



^ The Keldysh-Schwinger formalism has also been used in e.g. Ref. fl9\ in a non- 
thermal context. 
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p p 



\p-p'+pf 




Fig. 4. Graphical illustration of the r.h.s. of Eq. (16) 
p' 



P 



_ P 






Fig. 5. Amplitude A of Fig. |3] viewed as a dipole scattering amplitude. 

A can be viewed as a dipole scattering amplitude, as shown in Fig. |5j with the 
quark and the antiquark momenta, and respectively, given by 



P = {Po,Px,Py,Pz) = (p,0,0,p) 



P 



P±,P 



(17) 



The arrows on the quark lines in Fig. |5] refer to the flow of baryon number 



while the momenta are all directed to the right. The discontinuity in Eq. (16) 



is not just the imaginary part of the dipole scattering so that, in general, one 
cannot expect a simple relationship between the dipole scattering in Fig. [5] 
and the particular discontinuity illustrated in Fig. |4} 



3 Leading-order calculation 



In this section, we review the correspondence between the evolution of p±_- 
broadening amplitudes and dipole cross sections at leading order. Our aim is 
also to introduce the basic formalism and conventions that will be used for all 
calculations throughout this paper. 

Let us explain the method to evaluate a given lightcone perturbation theory 
graph such as the one shown in Fig. [6j There are essentially four elements: 
the product of the energy denominators D, the vertex factors, the gluon po- 
larization tensors, and the overall integration. 

In order to compute the energy denominators, we first assign a lightcone time 
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Fig. 6. Graph contributing to the leading-order p^-broadening amplitude. The 
freed gluon is emitted in the initial state in the amplitude and in the final state in 
the complex-conjugate amplitude. At the level of the nucleus, we chose one elastic 
interaction (momentum transfer li) and one inelastic one (momentum transfer 12). 
The exchanged gluons are drawn as dotted lines in order to simplify the graph. 

t to each vertex. We then evaluate D using the following formula: 

O = J Tl (f^ g*(^outgoing~^incoming)t±£i (18) 

where the product goes over all vertices and the integration takes into account 
the ordering of the different times, e is a small positive parameter (with di- 
mensions of an energy) that acts as an adiabatic regulator: It is set to zero at 
the end of the calculation. The sign in front of it is chosen in order to insure 
the convergence of the corresponding time integral at (plus or minus) infinity, 
depending whether the considered vertex is in the initial or in the final state. 
The energy of a state is the sum of the energies of all the particles in the state. 
The energy E of a particle of 4-momentum p = (j9+ ,p-,p) reads 

E = ^. (19) 

The quark-gluon vertices are evaluated in the eikonal approximation through- 
out this paper. For incoming/outgoing fermions of helicities a, a', colors a, b, 
momentum p, and a gluon which carries the Lorentz index n and the color C, 
the vertex is represented by the factor 

±2tgtZp,S''"\ (20) 



where the sign is + for a quark and — for an antiquark. The triple-gluon 
vertices which couple the gluons in the wave function with the nucleus will 
also have to be evaluated in the eikonal approximation. When the fast gluon 
line carries the Lorentz indices a, /3, the colors A, B, and the momentum p, 
the vertex factor reads 

'^gfABCP^Qafi- (21) 



9 



The triple-gluon vertex must however be computed exactly when it appears 
in the quantum evolution of the wave functions. For incoming momenta pi, 
P2 and ps ordered counterclockwise with respective Lorentz indices a, /3 and 
7 and color indices A, B and C, the corresponding vertex factor reads 

gfABC [gaii.P'i - Pl)l3 + gpa{Pl - ^2)7 + g^p{P2 - P'i)o\ {2'K f5^{pi +P2+ Pa)- 

(22) 

At next-to-leading order accuracy, we will never need to consider four-gluon 
vertices. 

In addition, each leg of vertex through which the momentum p flows has a 
factor 

The components of the polarization tensor for the gluon are of 3 different 
types. For a gluon of momentum p = (p+,j9_, j9i,p2), the latter read 

P Pi 

d__{p) = , d_i{p) = — , dij{p) = -gij = 6,j {i = 1,2). (24) 
P+ P+ 

An integration over all + and _L components of the loop momenta may then 
be performed, although, as we will discover, a momentum-by-momentum iden- 
tification holds for some classes of graphs. The integration over the on-shell 
momentum k is performed as 

(27r)3 • 



We are now ready to address the calculation of the diagram depicted in Fig. |6j 
We interpret it as the representation of a T-product, like all other diagrams 



in this work. Using Eq. (18), we find that the energy denominators read 



D=^, (26) 
where E = {li — p)^/(2A;+), E = {I2 +p)'^/{2k^), while the polarizations read 

Y^d.Si -P) ■ d^, ■ d,4-h-p) = + (27) 

Including the vertex factors, we get, averaging over the quark color and helic- 
ity: 

^^NlC^l^S^^P^^ , (28) 

2 {p-mp+kyk^ 

where Cp = {N^ — l)/(2Ai"c). The transverse momentum conservation is en- 
forced with the help of the additional factor 6'^{k + p — h). 
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Fig. 7. Graph contributing to the lowest-order quantum correction to the dipole 
forward scattering amphtude. Here there are 2 dipoles interacting with the target 
nucleus. This graph corresponds to the p^-broadening graph in Fig. [6j 

From the T-product, we may compute the contribution of the particular graph 
represented in Fig. [6] to the distribution of transverse momenta of the quark. 
To this aim, we need to integrate over the 3-momentum of the gluon k and to 
convolute with appropriate unintegrated gluon distributions which appear in 
the description of the scattering with the nucleus. The factor which comes with 
each two-gluon exchange with a particular nucleon in the nucleus involving 
the transfer of the momentum / reads 

±/ff^pL (29) 

where the sign is +1 if the nucleon scatters elastically and —1 if it scatters 
inelastically. When the exchanged gluons attach to the same (quark or gluon) 
line, a combinatorial factor | must be taken into account. Another overall 
combinatorial factor takes care of the time ordering of the scatterings off the 
different nucleons. (For the particular graph in Fig. [6| none of these factors is 
needed). All in all, we get: 



(Pp 



ttsNc /■+°° dk. 



graph in Fig. [6] 



iV2 - 1 7o k+ 

d^hd%{p-h)-{p + k) 



asXg{x,ll)pL asXg{x,l^)pL . (30) 

A priori, we need to sum over any combination of elastic and inelastic ex- 
changes. Actually, they would all look like the ones shown in Fig. [6] and bring 
factors which are identical in the p^-broadening case and in the dipole case. 
Hence it is enough to consider a small number of scatterings: We need only 
make sure that the case we single out be sufficiently general. 

We can compute in the same way the dipole graph in Fig. [7j The only differ- 
ences with respect to the previous case are: (i) the energy denominators come 
with a minus sign, (ii) the antiquark-gluon vertices have opposite sign with 
respect to the quark-gluon vertices and (Hi) the scatterings with the nucleus 
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-OC' 




c d 

Fig. 8. Full set of the lowest-order corrections to the p^-broadening amplitude. 
Graph c is the one already represented in Fig. [6j and graph d has the same expression 
since it is related by a simple symmetry. 



all come with an overall "+" sign in Eq. (29). All the extra minus signs com 



pensate, and we find an exact momentum-by-momentum identity between the 
expressions for the graphs of Fig. [6] and [7} 

We see that we may adopt the following rule as for the relative signs between 
Px-broadening amplitudes and dipole amplitudes: Since the extra minus sign 
which comes from the gluons attaching to the antiquarks is compensated by 
the change of sign of the discontinuity of the scattering with the nucleus when 
one goes from p± broadening to dipoles, it is enough to take into account only 
the extra minus signs stemming from the couplings of the gluons in the wave 
functions to the antiquark. 

We have seen in detail how the particular graphs in Fig. [6] and [7] are related. We 
need to extend the discussion to all possible graphs. First, we note that the p±- 
broadening graphs in which the gluon is not produced in the final state trivially 
correspond to dipole graphs for which the quark or the antiquark is a spectator. 
So we shall not discuss this case. As for the form of the interaction with the 
target, we choose to consider one elastic and one inelastic scattering. The gives 
the most general momentum fiow, and since all vertices are eikonal, it would 
be straightforward to extend the discussion to any form of the scattering. 

In the case of p_L-broadening, there are four different nontrivial graphs, depend- 
ing on whether the gluon is emitted after or before the time of the interaction, 
in the amphtude and in the complex conjugate amplitude. They are depicted 
in Fig. |8} The full set of dipole graphs at leading order is represented in Fig. 19} 



We have computed graph c (Fig. [6] and [s]) and shown that its expression is the 
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same as the expression for the dipole graph 7 (Fig. [7] and [9]). Graph d is similar 
to c and it is straightforward to check that its expression is the same as the 
expression for the dipole graph in which the gluon attaches to the antiquark 
in the amphtude and to the quark in the complex conjugate amplitude. In 
other words, 

c = 7 and d = S. (31) 



We turn to graph a in Fig. |8j The scattering occurs after the emission of 
the gluon both in the amplitude and in the complex conjugate amplitude. 
The gluons exchanged with the nucleus may be chosen to scatter with the 
quark line only: It is easy to check that the graphs in which the fast gluon 
scatters cancel between each other. The polarization factor for the gluon is 

then d (— ^2 —p)- We keep two scatterings with the quark, one elastic and 

one inelastic. Leaving out the factors describing the interaction with the target 
and the integration over k+, the expression for the graph reads (compare to 



Eq. (|28])) 

a = -2/C| (32) 

The corresponding dipole graphs would be a, a', /3 and (3' in Fig. |9] We see that 
these graphs are purely virtual: Therefore, for the correspondence to work, it 
is crucial that all interactions with the nucleus happen with the quarks and 
not with the gluon in the p_L -broadening case, which is effectively verified since 
scatterings with the gluons cancel among themselves. 

For this process, there are also graphs with instantaneous gluon exchanges 
(graphs a' and /3'). We take them into account by modifying the d compo- 
nent of the polarization tensor in the corresponding causal graphs {a and /3 
resp.) as follows [14]: 

d_4k) ^ d_4k) (1 - ) = d__{k) X F , (33) 

where -Egiuon is the energy carried by the gluon, and S is the factor in the 
energy denominators associated to the emission of this gluon. The second 
equality defines the F-factor. Since the graphs containing gluons which are 
exchanged instantaneously may be effectively taken into account through the 
above modification in the expressions of the causal graphs, we will not draw 
them systematically in the following sections. We shall just perform the sub- 



stitution (33) whenever the kinematics allow for instantaneous exchanges. 



Let us evaluate the dipole graph a in Fig. [9j The energy denominators read 

D = 7^ r , (34) 

2ie{E-ie) ' ^ ' 
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where E = (p + /2)^/(2/c+). The divergence encoded in the 1/e singularity may 
be traced to the absence of time scale for the exchange of the gluon between 
the quark and the antiquark, making it equally likely for all times between 
t = — oo and t = 0. The expression for this graph reads, before integration 
over the momenta, 

a = -i9n^^-^. (35) 

2^ ^ kl 2ieE-ie ^ ^ 

and hence this graph gives a divergent contribution. Adding graph a' turns out 
to cancel the divergence (a fact which has already been noticed by Kovchegov 
\23\). To take to latter into account, we multiply a by the factor F defined in 



Eq. (33), the energy of the gluon being E, and J2 = E — ie: 



E — is ie , ^ 

F„ = l-^ = j. (36) 

The sum of the graphs a and a' reads 

a + a' = -g'C%—^^, (37) 

which is now finite. The sum of the graphs /3 + /3' has exactly the same ex- 



pression. Comparing to Eq. (32), we see that the following identity holds: 



a = a + a' + P + P'. (38) 

The purely final state graph h corresponds, as we may check by explicit cal- 
culation, to the complex conjugate graphs: 

6 = a + a' + /3 + (39) 

This completes the proof that at leading logarithmic order, the evolution kernel 
of the px-broadening amplitude is identical to the evolution kernel of the dipole 
^-matrix. 



4 Proof of the equivalence at next-to-leading order 



We consider systematically all relevant graphs in the view of proving the equiv- 
alence of the evolution kernels at next-to-leading order, namely, that the iden- 
tity ([T| is preserved by the energy evolution at that accuracy. Therefore, we 
keep treating the quark-gluon vertices in the eikonal approximation, but the 
triple gluon vertices which contribute to the evolution kernel are computed 
exactly. 
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Fig. 9. Set of leading-order virtual (left) and real (right) corrections to the dipole 
amplitude. In the graphs a and /3, the exchanged gluon is causal, whereas in a' and 
/3', it is exchanged instantaneously. Graph 7 is equivalent to the graph in Fig. [7| 
and 5 has the same expression. 

We classify the graphs according to the number of quark-gluon vertices. There 
are two, three of four such vertices. In the first case, the diagrams contain a 
gluon loop, but it turns out that we do not need to evaluate it explicitly to 
prove the equivalence. In the second case instead, the three-gluon vertex has 
to be computed and the different gluon polarizations fully discussed. Further- 
more, we find that a momentum-by-momentum identification no longer holds, 
and the identification results from an analytical continuation of the integra- 
tion over the longitudinal component of the momenta. Finally, in the third 
case, the classes of diagrams that have to be considered together for the iden- 
tification to happen are larger. Moreover, while it is enough to analyze the 
energy denominators in order to prove the equivalence, which makes each cal- 
culation quite easy, the graphs are numerous, and analytical continuation is 
also needed for some diagrams. 



4-1 Two quark-gluon vertices 



We will first analyze the p^-broadening graphs, then the dipole graphs, and 
the comparison will be drawn in the last subsection. 



4.1.1 p±-broadening graphs 

With two quark-gluon vertices exactly, the next-to-leading order graphs con- 
tain a gluon loop. One such graph contributing to p± broadening is represented 
in Fig. [lOj The other graphs only differ from that one by the chronology of 
the interactions. 
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Fig. 10. One particular graph contributing to p_|_-broadening at next-to-leading 
order and exhibiting two quark-gluon vertices exactly, as far as quantum evolution 
is concerned, with one elastic and one inelastic interaction with the target. The 
energies and momenta of the particles are labeled explicitly. 

Hence we group the graphs according to the time at which the quark-gluon 
and gluon-gluon branchings occur, in the amphtude and in the complex conju- 
gate amplitude, relatively to the time at which the interaction with the target 
nucleus occurs. Each such class of (nonvanishing) p^-broadening graphs corre- 
sponds to a particular set of dipole graphs, which we will discuss in Sec. [4.1.2 
after having reviewed all p^-broadening graphs. 



We anticipate the fact that the only nontrivial factors to compare between 
the contribution of a given set of p^-broadening graphs and of dipole graphs 
are the energy denominators: One may therefore avoid the discussion of the 
gluon loop. 

We shall name the p^-broadening graphs with the help of the ordering of their 
vertices with respect to the scattering. For example IFIF (Fig. 11) means one 
qg vertex at early times and one 3-gluon vertex at late times in the amplitude 
(left of the cut), and the same in the complex-conjugate amplitude (right of 
the cut). 



Both quark-gluon vertices in the initial state. Let us first consider the 
case in which the qg branchings are in the initial state both in the amplitude 
and in the complex conjugate amplitude (namely the qg vertices occur at 
negative times). The gg vertices may then occur either at positive or negative 
times. 

The first case we examine is when the gg branchings both occur in the final 
state. The full set of graphs is shown in Fig. ITTj We see that we can label 



all momenta exactly as in Fig. 10 since the topology of all the graphs is the 



same. Therefore the expressions for these graphs will differ only in the energy 



denominators. Let us evaluate the latter using Eq. (18). 



We label the time at the branchings as < for the qg and gg branchings in 
the amplitude respectively, and td < tc for the same branchings in the complex 
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IFIF 



IFFI 



IIFF 



Fig. 11. Complete set of graphs for which the qg vertices are in the initial state 
and the 3-gluon vertices in the final state. The interactions with the nucleus are 
understood since they are the same as in Fig. [To] for all three graphs. 

conjugate amplitude. The energy denominators corresponding to the graph in 



Fig. \W\ (or equivalently the leftmost graph in Fig. 11) are obtained from the 



evaluation of the following integral (see Eq. 18): 



-^IFIF 



r+co rO 

dta I dtf, / dt 

Jo J+oo Jo 



dtn 



exp [i [Eta + (^1 + E2-E) ih - Q - Et^ +e{ta-h-t, + t^) } . (40) 



We have defined (see Fig. 10) 



E 



2k. 



E 



2k~ 



El 



2A;n 



E2 



2h 



2+ 



(41) 



and from the conservation of the longitudinal momentum, /c2+ = k^ — ki^. 
(Of course, the "bar" does not mean complex conjugate in these equations). 
The result of the integrations over the different times reads 



D 



1 



IFIF ~ - / - - X 2 • 

EE ( El + E2 — E 1 



(42) 



We perform a similar calculation for the two other graphs in this group. The 



one in which the loop is in the amplitude (Fig. 11, middle) has the following 
energy denominators, after expansion for e — 0: 



D 



IFFI 



2i€ EE (Ei + E2-E) 2 

1 

+ 



EE ( El + E2 — E ) 



E^E [El + E2-E) EE^ (El + E2-E)\ 



■ (43) 



This time, there is a divergent term when e — )■ 0, which is however imaginary. 
The graph in which the loop is in the complex conjugate amplitude has a 
similar expression up to sign changes, in such a way that in the sum -Djppj + 
Djjpp, the divergent imaginary part cancels as well as the two last terms in 
Eq. (43 ) while the second term gains a factor 2. The sum of the 3 denominators 
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h-p 




IIFI 





IIIF 



Fig. 12. Graphs in which the two quark-gluon vertices are in the initial state, one 
3-gluon vertex is in the initial state in the amplitude, and the other one is in the 
final state, either in the amplitude or in the complex-conjugate amplitude. This set 
should be supplemented with the symmetric ones, obtained by taking the complex 
conjugate. 



then cancels, 



D 



IFIF 



+ ^IFFI + D 



IIFF 



(44) 



and therefore, the total contribution of these graphs is zerop] 



Let us turn to the case in which the leftmost gg branching occurs before the 
time of the scattering in the amplitude while the rightmost one occurs in the 
final state, either in the amplitude or in the complex conjugate amplitude. 
Then there are the four diagrams depicted in Fig. [12] The Lorentz structure is 
the same for the graphs IIFI and IIIF. We evaluate the energy denominators 



of the causal graphs (first row in Fig. 12 ) as 



D 



IIFI 



EE{Ei + E2)iE[ + E'2- E') 



(45) 



where 



Eo 



{p- h + kl 



2h 



2+ 



{ki-hf 
2h+ 



E' 



{p + h + hf 



2k 



E' 



2+ 



2k. 



(46) 



see Fig. 12 Computing in the same way -Djjjp, we find that -Djiip 



-D 



IIFI- 



Since all other factors are the same for these graphs, their sum cancels, inde- 
pendently for each polarization of the gluons. 



^ Such cancellations, which are related to probability conservation and occur when 
one sums over a sufficiently inclusive subset of graphs, appear in many different 
contexts, see e.g. Ref. [22]. 
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nil nil nil 

Fig. 13. Causal graphs in which all vertices are in the initial state. The scattering 
with the target occurs with the quark: Graphs in which one of the gluons scatters 
cancel among each other. 




Fig. 14. Instantaneous-exchange graphs corresponding to graph IIII in Fig. 13 



They are effectively taken into account in our calculation by multiplying IIII for the 
longitudinal polarization of the gluons which attach to the quarks by two F-factors 



(see Eq. (33)). 



When the polarization ( ) is chosen for the leftmost gluon in IIFI and IIIF, 

then each of these graphs also cancels with the same graphs where the latter 
gluon is exchanged instantaneously (graphs in the second row in Fig. 12). 
Indeed, according to the rule (33), the polarization tensor d for this gluon 



has to be multiplied by the factor 



-^IIFIJIIF — ^ 



E 



le 



E 



E' 



(47) 



The product of this factor with the corresponding energy denominators van- 
ishes since -DufijiIf finite for e — )■ 0. So there is no global contribution of 
the instantaneous-exchange graphs either. 

So far, we have not found any contribution to p^-broadening of the graphs that 
we have examined. We turn to the last class of graphs in which both three- 
gluon vertices are in the initial state, either in the amplitude or in the complex- 
conjugate amplitude. Again in this for the leading-order graph a in 
Fig. |8| all graphs in which one of the gluons in the wave functions scatters with 
the nucleus cancel among themselves. The complete set of graphs in which the 



gluons are causal is shown in Fig. 13 There are in addition graphs in which 



at least one gluon is replaced by a contact interaction, see Fig. 14 and 15 



Before we proceed to the computation of the energy denominators, we need 
to discuss the helicity structure of the gluons. The leg of the gluons which 
attached to the quark always carries the longitudinal polarization "— " since 
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X -K 



•a X 



Fig. 15. Instantaneous-exchange graphs corresponding to graphs IIII and IIII in 



Fig. 13 



the qg vertices are eikonal. The other leg of these gluons is either longitudinally, 
or transversely polarized. Actually, in these graphs in which no interaction 
occurs between the gluons in the wave function and the target, the leg which 
attaches to the loop of both of these gluons must have the same polarization. 
Hence both gluons come either with a d or with a polarization factor. 

In order to easily incorporate the instantaneous-exchange graphs using the 



modification (33) of d , we will distinguish the two possible polarizations: In 

the (— -L) case, it will be enough to consider the sum of the energy denomi- 
nators of the relevant graphs, while in the ( ) case, we will have to weight 



them by appropriate F-factors computed with the help of Eq. ( 33 ) 



We go back to the graphs of Fig. 13 The energy denominators read 



The sum of the denominators reads 



Dim + Am + Am - ^2(^^ + ^ E\Ei + E^) " ^^^^ 



The incorporation of the graphs of Fig. 14 and 15 require to multiply -Dug 
and D^YTi by 

E — is is is 

and Djjjj by Fjjjj x -Fjjjj. So because of the finiteness of the energy denomina- 
tors, once we add the instantaneous-exchange graphs, the contribution of the 
( ) polarization vanishes. 



(50) 



Both quark-gluon vertices in the final state. There are 3 causal graphs 
in which the gluon is emitted off the quark at late times both in the amplitude 
and in the complex conjugate amplitude. They are represented in Fig. 16 Since 
the kinematics is the same for all these graphs, it is again enough to address 
the energy denominators. 
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X 



FFFF FFFF FFFF 

Fig. 16. Set of the causal graphs in which all vertices are in the final state. One 
should also consider the graphs where some gluons are exchanged instantaneously: 



We do not show them since they look exactly like those of Fig. 14 and 15 
They read 

1 

D 
D 



^^^^ (^1+^2)2(^1 + ^2-^)' 



1 



FFFF 



2ie {El + E2-E) E^ {El + E2 - E) (51) 
1 1 



2 E^ {El + E2- Ey 



-^FFFF — -^FFFF- 



After some rearrangements, their sum reduces to 



1 2 

^FFFF + ^FFFF + ^FFFF = E^Ei + E^Y ^ E\Ei + E^Y ^^^^ 



namely it is identical to Eq. (49 ). In order to take into account the instantaneous- 
exchange graphs, we need to multiply the denominators by the respective 
factors 



-^FFFF — ( ) ' -^FFFF ~ ~^ £ "^FFFF ~ -^FFFF' (^3) 



Then we find that the sum of all instantaneous-exchange graphs and of the 

causal graphs with the ( ) polarization for the gluons that couple to the 

quark has the factor 



■^FFFF-^FFFF + "^FFFF-^FFFF + "^FFFF-^FFFF ~ ^i (^^) 



and hence vanishes. 
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FFFI 



FFIF 



FIFF 



Fig. 17. Causal graphs in which one qg vertex is in the initial state, and one in the 
final state. This list should be completed with the instantaneous-exchange graphs 
in correspondence with the two leftmost graphs, and with the complex-conjugate 
ones. 

One qg vertex in the initial state, one in the final state. The causal 
graphs are represented in Fig. 17 The energy denominators read 



D 



1 



1 



FFFI 



+ 



1 



2ie EE{E - El- E2) 2EE^{E - E^- E2 



+ 



+ 



for the leftmost graph 
D 



2E^E{E - El- E2) 2EE{E - El- E2Y 



1 



(55) 



FFIF 



E(^Ei + E2){E — El — E2 



(56) 



for the graph in which there are 2 gluons in the final state (graph in the middle 
in Fig. 17), and 

1 



D 



1 



1 



1 



^^^^ 2ie EE{E - Ei - E2) 2 EE^{E - Ei - E2) 

11 11 



(57) 



2E^E{E - El- E2) 2EE{E - El- E2Y 

for the graph in which the loop is on the right of the cut (rightmost graph in 
Fig. 17). The sum of these denominators reads 

-^FFFl + -^FFIF + -^FlFF ' 



1 



E^E{Ei + E2y' 



(5J 



The instantaneous-exchange graphs are taken into account with the help of 
the factors 

2ie _ El + E2 

-^FFFI — 1 + -^FFIF ~ ^ ' i^^J 

Again, the contribution of these graphs summed with the causal graphs which 
have the ( ) polarization vanishes, since 



■^FFFI-^FFFI ~^ -^FFIF-^FFIF ~^ -^pttttt — 0. 



'FIFF 



(60) 
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Fill 

Fig. 18. Causal graph in which one qg vertex is in the initial state, one in the final 
state, and the gluon loop is fully in the initial state. 




h — p 




h-p 



FIIF 



FFII 



Fig. 19. Causal graphs in which the gluon loop interacts. There are in addition 2 
instantaneous-exchange graphs (not shown). 



The graph in which the gluon loop is in the initial state (Fig. 18) has the 
following denominators: 



D 



1 



Fill 



EE^{Ei + E2)' 



The sum of the instantaneous-exchange graph with the 
the causal graph once again vanishes. 



(61) 

polarization in 



Finally, we address the graphs in Fig. 19 We find the following energy de- 
nominators: 



D 



FIIF 



-^FFII 



1 

EE{Ei + E2){E-E'{-E'^' 
1 



E{Ei + E2){E'i + E'^{E - E'{ - E'^ 



where 



E'i 



,„ {ki + h + hf ^„ (p + h + h) 



2k 



1+ 



, E2 



2k 



2+ 



Their sum reads 



(62) 



(63) 



FIIF ~ -^FFII 



EE{E, + E2)iE'{ + E!^) 



(64) 
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Fig. 20. Virtual graphs (One dipole exactly scatters). There would be 6 additional 
graphs in which one or both gluons which couple to the quark/antiquark are ex- 
changed instantaneously, ai and a2 only differ by the ordening of the qg and qg 
vertices. We do not display the complex-conjugate graphs. 

Again, the sum of the ( ) components and of the instantaneous-exchange 

graphs gives a null contribution, due to the finiteness of the denominators and 
to the fact that the F-factors are of order e. 



4- 1-2 Dipole graphs 



The full set of nontrivial causal dipole graphs is shown in Fig. 20|21||22 up 



to graphs deduced from the latter by obvious symmetries. (We do not discuss 
the graphs in which the gluons both couple to the same quark or antiquark, 
because the correspondence with p^-broadening graphs in which the evolution 
happens entirely either in the amplitude or in the complex-conjugate ampli- 
tude is then immediate). We may label the flow of the momenta through the 
graphs in such a way that it is the same as in the p^-broadening case. This 
justifies a posteriori that only the energy denominators need to be compared. 

The expression for the energy denominators are very easy to obtain for these 



graphs. The graphs a in Fig. 20 are virtual graphs, which contribute equally. 
The energy denominators read 

_ 1 1 2 1 

Da, + - - — ^^^^—^ - —^——^ - ———^ (65) 

where E, Ei E2 are defined in the previous section. To compute the sum of 

the ( ) polarization and of the instantaneous-exchange graphs (not drawn), 

it is enough to multiply the previous denominators by the factor 

E -ie\ f E - ?,ie\ , , , 

1 - (1 ^) ~ e\ (66) 

and thus the latter sum vanishes: The ( ) polarization of the gluons does 

effectively not contribute. 
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Fig. 21. Graphs where there are 2 dipoles at the time of the interaction. There 
would be 2 additional instantaneous-exchange graph, and also the graphs in which 
the gluon couples to the quark in the amplitude and to the antiquark in the com- 
plex-conjugate amplitude (namely the complex-conjugate graphs). The arrows on 
the gluons indicate the direction of the flow of the transverse momentum. 




7 

Fig. 22. Graph where there are 3 dipoles at the time of the interaction. In the 



same manner as for the diagrams of Fig. 20 there would be 3 extra graphs with 



contact interactions instead of the external causal gluon propagators. There are 
also the graphs in which the couplings of the gluons to the quark and antiquark are 
interchanged between the amplitude and the complex-conjugate amplitude. 



Similarly, the energy denominators for the graphs of Fig. 21 read 



E^E{Ei + E2f' EE^{E^ + E2Y' ^^^^ 



Last, the graph in Fig. 22 gives 



EE{E^ + E^){E'l + E'^' 



(6J 



The F factors for the two last classes of graphs are proportional to e: Hence the 

( ) polarization components compensate with the instantaneous-exchange 

graphs also in the two previous cases. 
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4.1.3 Proof of the correspondence 



The identification between the relevant graphs is now quite straightforward. 
The nonzero graphs in the case of p^-broadening are those of Fig. 13} 16, 17 



18 , 19 with the 



polarization for the gluons which hook to the quark. The 



associated denominators are given in Eq. (49), (52), (58), (61), (64) respectively. 



We naturally identify the dipole graphs ai, ^2 of Fig. 20 with the p^-broadening 



graphs where all vertices are in the initial state of Fig. [13} since they can be 
topologically related by bending over the quark line to form the antiquark of 
the dipole. It is easy to check the formal identity of the expressions for these 
graphs: The polarization and vertex factors are the same in the p_L-broadening 
case and in the dipole case, except for a minus sign difference which comes 
from the qg coupling becoming a qg coupling. We notice that we may also 
choose the interaction with the nucleus to be the same in both cases, and 
label the momenta in the same way as for p^-broadening. We see that the 
following identity between the energy denominators holds: 



D 



iiii 



D 



IIII 



D 



mil 



-Re {D^, + 



(69) 



The minus sign in the r.h.s. is explained by the coupling to the antiquark. 
The only remaining mismatch is an imaginary term which cancels when the 
complex conjugate graphs are taken into account. 



Bending over the quark on the right of the cut in the graphs of Fig. 16, we see 



that we get the symmetric of the dipole graphs of Fig. 20 , namely the complex 
conjugate graphs (not represented). The formal identity between the energy 
denominators reads 



D 



FFFF 



~l~ ^FFFF ~^ ^FFFF — (-^cti ~^ -^05 



(70) 



In the same manner, we see that the graphs of Fig. 17 are topologically related 



to f3i in Fig. 21, the one in Fig. is related to (32, and the ones in Fig. 19 



look equivalent to 7 in Fig. [22} As a matter of fact, the following identities 
are verified: 



D 



FFFI 



+ Dppjp + D 



FIFF 

-^Fm 



-D 
-D 



(71) 



D 



FIIF 



+ D 



FFII 



Equations (69),(70),(71 ) prove the equivalence of p± broadening and dipoles 



for the considered class of diagrams. 



Note that the p^-broadening graphs which cancel among themselves would not 
have any dipole counterparts allowed by the kinematics. Indeed, for example 
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X 



Fig. 23. Dipole graph which would correspond to the (vanishing) p_|_-broadening 



graph of Fig. 12 , but which is not allowed by momentum conservation. 




Fig. 24. Typical momentum configuration around the 3-gluon vertex. The arrows 
show the direction of the transverse components of the momenta. The "+" compo- 
nents are always directed to the right. The gluons which carry the momenta ki and 
/cs may also be exchanged instantaneously in some graphs, in which case the factors 
El and which come with the polarization tensors of these gluons (see Tab. [T]) are 
replaced by Ei — Si and E^ — T,^, where Si and S3 are the energy denominators 
of the corresponding causal graphs. 



the graphs of Fig. 12 would correspond to a dipole graph such as the one 



shown in Fig. 23, which is forbidden by momentum conservation. 



4-2 Three quark-gluon vertices 



We turn to the case in which there are 3 quark-gluon vertices. The NLO 



graphs all have one 3-gluon vertex (like the one in Fig. 24), see e.g. Fig. 25 for 



a sample of the p_L -broadening graphs and Fig. [26] for the would-be equivalent 
dipole graphs. 

Let us analyze the Lorentz structure of the upper part of these graphs. The glu- 
ons entering the vertex may have "— " or "±" polarizations. It turns out that 
there are only four possible configurations: Either all gluons are transverse, or 
one of them (and only one) has a "— " polarization. All other possibilities lead 
to vanishinig contributions. The gluon which is polarized longitudinally may 
be exchanged instantaneously if the kinematical configuration of the diagram 
allows it; Such configurations are taken into account by a modification of the 



corresponding polarization tensor, see Eq. (33) 



It is useful to write down the factors specific to the different polarization config- 



urations for the typical graph shown in Fig. 24 All other configurations of the 



momenta may be deduced from the one illustrated in Fig. p4^ by appropriate 
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Vertex and polarization factors 



Vi 
V2 



-Ei{k2 + h) + {k'^-k^) 
m^{k2-k',){k,.k,) 

E^{ki + k2) + {ki-A) 



{k2 ■ k3){k^ ■ k'^) - (ki ■ k2)ik3 ■ k'^) 



Table 1 

Factors V associated to the 3-gluon configuration shown in Fig. [24] for the different 
polarization configurations of the gluons at the 3-gluon vertex (see the main text). 
The transverse momenta are always incoming, while all "+" components of the 
momenta are flowing from the left to the right. If k2 = k2, then V± = and 
V2 = 2E2{ki - ks)+h ■ h/{ki+k2+k3+). 



changes in the signs. Note that we will keep the transverse momenta always 
incoming to the 3-gluon vertex in order to simplify the comparison of the 
graphs, at the expense of having sometimes the "+" and "_L" components of 
the momenta flowing in opposite directions (in which case d-±{k) = —k/k+). 



We contract the 3-gluon vertex with the gluon polarization tensors assuming 
that the 3 external legs of the gluons represented by wavy lines couple eikon- 
ally, thus always have the "— " polarization. The gluons which are exchanged 
with the target (represented here by the dotted vertical line in Fig. 24) also 
couple eikonally to the leftmost gluon. We temporarily leave out the color 
factors, the coupling constants, and the momentum factors (23) that come 
with each leg of vertices since the latter may be put back in the end and do 
not depend on the polarization. There are four possible polarization configu- 
rations that we need to distinguish. Either one gluon (and only one) has the 
"— " polarization at the 3-gluon vertex, or all have the "_L" polarization. To 
each of these configurations corresponds a specific factor V. 



We present the factors V for the four possible polarization configurations in 
Tab. [Tj The index of V corresponds to the label of the gluon which enters 
the 3-gluon vertex with a minus polarization, or is "_L" when all gluons are 
transversely polarized. 



Coming back to the full diagrams, we have to distinguish the different possible 
time orderings of the vertices. Some orderings lead to real dipole graphs, that 
is, where there are two dipoles in the wave function at the time of the inter- 
action, while other orderings lead to virtual dipole graphs, which contribute 
to the renormalization of the wave function of the initial dipole. 
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A 



B 



Fig. 25. p_L-broadeiiing graphs in the case in which two quark-gluon vertices are 
in the initial state and one is in the final state. The transverse momenta are di- 
rected according to the arrows, while the "+" momenta are oriented to the right. 
We do not show the graph which is similar to A but with gluon 3 exchanged 
instantaneously, and -B^, similar to B but with gluon 1 exchanged instantaneously. 
Momentum conservation implies that k'2 = k2 + h + h- 




a 





4 


























P 




'ii! it J 








-k—k 





/3 




7 



Fig. 26. Dipole graphs corresponding to the p^-broadening graphs in Fig. 25 To this 
list one should add the graphs in which gluons 1 or 3 are exchanged instantaneously, 
and distinguish two relative orderings of the leftmost qg and qg vertices in 7, which 



we should draw as two different graphs 7^"^^ and 



.(2) 



4-2.1 Initial state/initial and final state 



Let us start with graphs related to real dipole graphs, and first with the 
simplest set of them, characterized by two of the qg vertex and the 3-gluon 
vertex being in the initial state, while the remaining qg vertex is in the final 



state (see Fig. 25). On the dipole side, the corresponding process involves two 



dipoles at the time of the interaction (Fig. 26). 



The most general case for the interaction with the nucleus is captured by two 



scatterings such as in Fig. 25 Higher-order scatterings add simple combina- 
torial factors. 



The calculation of the different contributions using the previously edicted 
rules (Eq. 



|l8j) for the energy denominators, Eqs. (20) and (21 ) for the eikonal 
vertices, Eq. 23 for additional factors associated to all vertices, and Tab. [T] for 
the 3-gluon vertex dressed with the polarization tensors) is straightforward. 



29 



We define the energy of tlie gluons as 



El 



2k 



1+ 



2k^4 



E' 



h.12 

tv2 



2k 



2+ 



(fcs + ri + fa)- 

2k2^ 



2h 



3+ 



(72) 



Tlie calculations corresponding to the A, a and B, /3 graphs (without includ- 
ing the lower part of the graphs representing the coupling to the target) are 
summarized in Tab. [2] and 3j We distinguish the four different polarization 
configurations (as in Tab. [T) for each graph. The two columns single out the 
energy denominators (which are not polarization dependent) and the 3-gluon 
vertex and polarization factors read off Tab. [T] with the relevant kinematics. 
In addition, there is an overall factor, which is identical for all graphs: After 
averaging over the color and helicity of the initial quarks, it reads 



{ki+k2+k^+Y' 

The 2/(A;i+A;2+/i;3+) factor in Tab. [T]has been incorporated in Eq. (73). 



(73) 



In order to compute the contribution of these graphs to p_L -broadening, we 



need to convolute the content of Tab. 2 and 3 multiplied by Eq. (73) with 



some weight which encodes the model for the interaction with the target and 
which can be computed from the rules listed in Sec. |3| 



dN 



contribution of a 
particular grapli 
■ + 00 

X 



dHi dH2 



-pL 



m - 1 



pL 



X 



r+oo 

/ dki^dk2+dk3^6{^^+" momentum) x [overall factor] 
Jo 

d^ki d^k2 d^k^ ^.i 



(27r)2 (27r)2 (27r) 



-5 (fci + /i;2 + fcs) 5 (global "_L" momentum) 
X [right column of the corresponding table] (74) 



where, for example in the case of graph A, the graph-dependent momentum 
conservation factors read 



5("+" momentum) = (5(fci+ + A;2h 

^^(global "_L" momentum) = b'^{k2 +p) 
[overall factor] = Eq. (73) 



-k. 



(75) 



and the "corresponding table" is Tab. [2] 

From the conservation of 3-momentum and from the positivity of the 
components, we note that the A and a graphs have the same integration 
range of the "+" components of the momenta, namely < < /C3+, while 
in the case of B and /3, k\j^ > k^^. As for the 7 graph, ki^ and k^^ must be 
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A, 



02 



an 



Energy denom. 



EiE'^E2,{E'i_+E2) 




E2E'^E-i(Ei+E2) 



3-gluon vertex xpola. 
{k'^-k^)Ei{2k^-ki)+ 
ih ■ h)ik2 ■ k'^) [-{ki + k^)+/{2{k^ - k^)+)] 


\k-h){ki-l)-{ki-k2){h-l) 



{k'^-k3){-E2)i2k3-ki) + 

{ki ■ ■ k'^) [-{ki + A;3)+/(2(A:3 - k{)+)] 



5 \{k2-k:i){ki-l)-{ki-k2){hi-l) 



Table 2 

Expressions for the graphs A and a in Fig. 25 and 26 , k'2 = k2 + l (where / = ^1 + ^2) 
and according to momentum conservation, k2+ = ^3+ — The "*" superscript on 
the graph labels means that an instantaneous exchange graph has been added. The 
factor 2/(fci+A;2+^3+) present in Tab. [T] was left out and integrated to the overall 
factor (|73|). The contribution of A to the p_L -broadening amplitude is obtained from 



the convolution given in Eq. (74). 



Energy denom. 



Bl 
B2 
B3 
B± 



ExE'^E3{E2+Ez) 



PI 

/3± 



1 

EiE2E'^{E2+Es,) 



3-gluon vertex xpola. 





(fci • k^){k2 ■ k'^) [-(fci + A:3)+/(2(A:i - ^^3)+)] 
{ki ■ k'2)Esi2ki - k3)+ 
-\ \{k2-k3){ki-l) + {ki-k2){h-l) 







{ki ■ k^){k2 ■ k'^) {{ki + A:3)+/(2(A:i - A;3)+)] 
{ki ■ k'2)E2{2ki - k3)+ 
\\{k2-k){ki-l)-iki-k2){k-l) 

Table 3 

The same as Tab. [2] but for graphs B and /3 in Fig 25 and 26, for which momentum 



conservation reads h 



'2+ 



k 



1+ 



k. 



both integrated from to +00. We may subtract the dipole graphs from the 
production graphs whenever the kinematics are identical. We find that only 
the "2" and "_L" components are nonzero in the differences A — a and B — (5. 
They are listed in Tab. |4] together with the nonzero components of 7. We have 
specified the integration range and singled out the factors which depend on 
the "_L" and "+" components of the momenta. Other factors are common to 
all graphs (see the caption of Tab. |4]). As for 7, we were careful to distinguish 
two possible orderings of the qg and qg vertices 7'^^'' and 7*^^^ before adding up 



their contributions, see the caption of Fig. 26 
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A2 — a2 
A± — a± 



B2 
B± 



P2 
/3± 



72 



7± 



Integration range 



< ^1+ < k?,+ 



'3+ 



< < +00 



Momentum-dependent factors 



(A.l-fc3)(fc2-fc^) 

ftj^/vj/tj '(■3 
-(fc2-fc3)(fcl-r) + (fcl-fc2)(fc3-r) 

t2 t2 p2 u2 
"■1^2 "'2 "■3 



(fci-fc3)(fc2-fcg 

^2 u2 p2 u2 
X 2 2 3 

-(fc2-fc3)(fcrO + (fcl-fc2)(fc3-0 
J:2 P p2 t2 

(k^_-k3){k2-k'^) 
t2 t2 p2 %2 

-(k2-k-i){kld) + {ki-k2){ki-l) 
t2 u2 p2 12 



fcl++fc3+ 
fcl+{fc3+-fcl + )fc3+ 



1 



-^34 



fci++fc3+ 
fcl+(fc3+-fel + )fc3+ 

1 

^1 + ^3+ 

fc3+-fcl + 
fcl+(fcl + +fc3+)fc3+ 

1 

fcl + fc3+ 



Table 4 

Differences between the contributions of graphs A,a and f3 respectively, and 
nonzero contributions of graph 7. We only keep the momentum-dependent factors, 
which we split between a factor depending only on the transverse components, and a 
factor depending only on the longitudinal ones. (The overall constant, which gathers 
the coupling and color factors, reads 647:'^a'^N^CF)- 



We are now in position to prove the equivalence between the production and 
the dipole processes. We need to compare the results of the integration of the 
expressions for the different graphs over the variable. We see in Tab. |4] 
that A — a and B — (5 have exactly the same expressions: The difference is 
in the integration ranges in which turn out to be complementary on the 
positive real axis [0, +00 [. 



Let us show in detail that the analytical continuation of 72 to the negative /ci+ 
region has exactly the same expression as 742 — 0:2 and B2 — P2, once the change 
of variable — t- — A;i+ is done in order to go back to positive + factors. To 
this aim, it is enough to focus on the "+" components in Tab. |4] (see the 
rightmost column) since the labeling conventions have been chosen such that 
the other factors remain identical for all graphs. We define x = ki^/k^^ in 
order to simplify the equations and write the integral associated to 72 as 



72 ^ 



X 



dx— r 

A Xil + X) 



(76) 



We go to the complex plane of the x-variable and sum the integrand over the 
double contour shown in Fig. |27| The integrals over the semi-circles /a, /a, le 



and their symmetries give a finite contribution, but their sum vanishes. From 
the Cauchy theorem, the remaining terms are then related by the equation 



2(/, 



0. 



(77) 
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Fig. 27. Complex plane of the x- variable in the integrand of Eq. ( 76 ) 
The first term J+ is proportional to 72. The two other terms read 



^ ^ 1 — X 

ax— r 

-i+e a;(l + X) 

-1-= 1-x 

-A Xl 1 + X) 



A 



ax^-. ; — > —{^2 



A 



x{l — X) 
l + x 

l+e X[l — X) 



[Bo 



(^2) 
P2) 



After the limit e, A — )■ and A — )■ +00 is taken, Eq. (77) shows the equiva- 



lence between the p^-broadening diagrams and the dipole diagrams. The same 
argument applies to the "_L" component and leads to the same conclusion. 



Note using the double contour in Fig. [27] is equivalent to taking a principal- 
part prescription for the pole at x = —1. Originally, this pole was shifted off 
the real axis by the regularization chosen for the energy denominators. But 
since the imaginary terms anyway cancel after one has included all relevant 
graphs, these two choices are equivalent. 

Of course, in the present case, the integrals /+, /< and /> are simple enough 
that we can also afford to check directly the identity 



k3+-£ 



dki+{A2,± - a2,±) + / dki+{B2,± - /32,j 



dki+-f2,±, (79) 



which holds once appropriate cutoffs have been introduced to regularize the 
infrared and ultraviolet divergences at ki+ — )■ and — )■ +00 respectively. 
But the analytical continuation argument presented above appears more gen- 
eral. 
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X— XX— X x-x 



Fig. 28. px-broadening graph for the configuration in which one quark-gluon cou- 
pling and the 3-gluon vertex are in the initial state, while the remaining two qg 
vertices are in the final state (left) and the equivalent dipole graph (right). 

Having performed these detailed checks in this simple case, we will be able to 
compare the graphs in all other cases by going through the following steps: 



(1) We write the expression for all graphs for all possible polarizations of the 
gluons, with conventions for the momenta that are such that the same 
transverse structure is kept for all graphs (This step is summarized in 
Tab. [2] and |3] in the case just discussed). 

(2) We choose 2 variables among the 3 "+" components ki^, ^2+, ^3+ and 
use momentum conservation in order to express one of them (/c2+ here) 
in terms of the other two. We fix one of the remaining variables (^3+ 
here). We group all graphs according to the range of integration of the 
leftover variable (fci+). In each such group, we add up the expressions 
for the px-broadening graphs and subtract the ones for the dipole graphs 
(Tab.g. 

(3) In the class of graphs for which the range of integration is unrestricted 
(graph 7 here, see Tab. |4]), we invert the sign of the variable on which we 
integrate (A;i+ here), and add an overall "— " sign to take into account 
the change of sign in the integration element. 

(4) The equivalence between p^-broadening and dipole amplitudes is true 
if the expressions obtained after this transformation are the same as the 



expressions found for all other classes of graphs (Eq. ( 79 ) in our example 



This procedure is a practical formulation of the analytical continuation method 
detailed above, which enables one to check at first sight the equivalence. 



4-2.2 Initial state/final state 

We now turn to the case in which one quark-gluon vertex in the amplitude is 
in the initial state, while the two other quark-gluon couplings on the right of 
the cut are in the final state. 

In the simplest case, the 3-gluon vertex is also in the initial state. This is 
shown in Fig. |28j This case is actually straightforward. On one hand, the 
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Fig. 29. p_L -broadening graphs for the configuration in which ah interactions but one 
quark-gluon couphng occur in the final state (The instantaneous-exchange graphs 
are not shown). 
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Fig. 30. Dipole graphs corresponding to the p^-broadening graphs of Fig. 29 (The 
instantaneous-exchange graphs are not shown). 



energy denominators, for px-broadening as well as for dipoles, read 



D 



1 



E'^Es{E, + E2)iE'{ + E'^y 
On the other hand, all other terms are identical for both graphs. 



^0) 



The graphs in which the 3-gluon vertex is in the final state are shown in Fig. 29 



(and Fig. 30 for the corresponding dipole graphs). We may group the graphs 



A, C and a in Fig. 29 and 30 since they share the same kinematics, while 
we treat separately the graphs B and (3. We do not present the evaluation of 
all diagrams separately for this case (step 1 in the procedure outlined above) 
since their expression is very similar to the previous case. We find that the 
expressions for A + C — a and for B are identical once we express everything 
in terms of the variables and k^^. 



The overall factor is again given by Eq. (73). The results for the nonzero 
components of A + C — a, B and (3 are shown in Tab. [5j 

It is clear that the analytical continuation of the components of the graph /3 in 
the variable ^2+ is equal to A + C — a for k2+ > ^3+ and to B for k2+ < ^3+ as 
far as the real part is concerned. The imaginary parts cancel with the complex 
conjugate graphs that also have to be taken into account. 
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Integration range 



Momentum-dependent factors 



pure "_L'' 



■'+" and "_L" 



^3+ < ^2+ < +00, 
< A;2+ < ^3+ 



(fci-fc3)(fc2-fc^) 



1 



A2 + C2-a2, B2 
A^ + Cl-al Bs 
A± + C±- a±, B± 

Table 5 

Sums and differences of the graphs of Fig. [29] and 30 grouped according to kinemat- 
ics. The overall constant is QAtt^o^N^Cf. 



< /C2+ < +00 



{k2-k3)(ki-l)-{ki-k2)(k3-l) 

J.2l.2r,/2 
''■l'''2'''2 

(fcl-fc3)(fc2-fc^) 

t2 t2 y2 
'^l'^2'^2 

(fcl-fc^) 

?2p2 
ft J Kg 

{k2-k3){ki-l)-iki-k2){k3-l) 
i.2l.2i,/2 
/(,-^«,2''-2 



fc2/fc2+_fc2/A;3+ {k2+-k3+)k2+k^^ 

1 2fc2+-fc3+ 

fc2/fc2+-fc|/A;3+ ik2+-k3+)k2+kl^ 

1 1 

kl/k2+-k'l/k3+ ik2-k3) + k'l^ 



fc2++2fc34 



k^/k2++k1/k3+ (fc2++fc3+)fc2+fcf!^ 
1 2fc2++fc3+ 

k^/k2++kl/k3+ ik2++k3+)k2+kj_^_ 

1 1 

kyk2++kl/k3+ (fc2+fc3) + fc^+ 




Fig. 31. Set of the -broadening graphs for which all interactions occur in the 
initial state but which cancel. We have not drawn explicitly the gluons exchanged 
with the target. 






B 



C 



Fig. 32. The same as Fig. [31] but for the non-cancelling graphs. The gluons ex- 
changed with the target (not drawn) all couple to the quark line, and thus are not 
relevant to the discussion of the equivalence. 

4-2.3 Initial state/initial state 

We turn to the graphs in which all vertices are in the initial state. First, 



there is the trivial case of the p^-broadening graphs in Fig. 31 which cancel 
in exactly the same way e.g. as the graphs of Fig. [12] cancelled. There is no 
kinematically-allowed equivalent dipole graph. 



We turn to the more interesting graphs in Fig. 32 All interactions between the 
nucleus and one of the gluons cancel, hence the exchanged gluons necessarily 



hook to the quark line. The corresponding dipole graphs are virtual (Fig. 33). 



The evaluation of the p^-broadening graphs of Fig. 32 is straightforward. The 
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5.^ 



C3 



Integration range 



< A;i+ < k 



'3+ 




< A;i+ < k3+ 



Energy denom. 



EiEs{Ei+E2r 



EiEj{E2+E3) 



EiE'l{Ei+E2) 



3-gluonxpola. 

(^2- 4)^1(2^3-^1)4 
(k^ . k3){-E2){h + h) 






{ki-k3)E3{ki + k3)^ 

{h ■ k2)E3{2ki - k3) 



ik2 ■ k3)Ei{2k3 - ki)+ 

{ki-k3)Ei{ki+k3) + 
{ki ■ k2)E3{k3 - 2ki) + 



Table 6 

Expressions of the graphs of Fig. 32 The momentum conservation reads k2+ = 
A;3+ — kij^ for graphs A and C, and k2+ = A;i+ — A;3-|- for graph B. These expressions 



must be multiphed by the overall factor (81) 





a 



7 



5^^^ (up) and 5^^^ (down) 



Fig. 33. Set of the dipole graphs corresponding to the p_L-broadening graphs of 
Fig. 32 5^^^ and (5^^^ differ by the ordering of the absorption times of gluons 2 and 



3. 7 actually represents 2 graphs, which differ by the order of the leftmost qg and 
qg vertices. 

results are presented in Tab. [6] 



The dipole graphs in Fig. 33 bring in a new difficulty with respect e.g. to the 
real graphs in Fig. [28|29p0j The energy denominators diverge, thus it is cru- 
cial to keep the regulator e throughout. The divergencies are not removed by 
instantaneous exchanges as was the case for other graphs, but they give imag- 
inary contributions which are cancelled when one adds the complex conjugate 
graphs. The results are shown in Tab. [7} 

The overall factor, including the averaging over the color and helicity of the 
quark, reads 

327r^a^iV,C| 

{k,^k2+k3+r ^ ^ 



which is identical to (73) in the large- iV^ limit [Cp ^ Nc/2) except for the 
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1*2 
73 



Integration range 



0<ki+ < /C3+ 



< A;i+ < +00 



Energy denominators 



1 



1 



4i£ (E2-3i£)(E3-ie){Ei+E2-2ie) 



1 



1 



4i£ (Ei-i£){E3-3i£){Ei+E2-2ie) 



its (Ei-ie){E3-ie){E2-2ie) 



3-gluon vertex xpola. 
ih ■ h){2i£ - E2){2h - ki). 
{h-h){-E2){ki + k3)+ 
{ki-k2){ie){k3-2ki)+ 
{k2-k3)i-Ei)i2ks-ki)+ 
{ki-ks){2ie-Ei){ki + k3)+ 
{ki-k2){2,ie){2ki-k-i)+ 
{k2-k3){is){2ks + ki)+ 
{ki ■ k3){3ie){k3 - ki)+ 
{h-k2){-ie){2ki + k3)+ 



(1)* 



^(1) 

"2 



(2)* 



4i£ {Ei-i£){E2+E3-2i£){E2-3i£) 



ks+ < ki+ < +00 



ik2 ■ k3)iie){2ks - ki)+ 
(ki ■ ks)E2{ki + ks)+ 
ih-k2)iE2-2ie){2ki-k3) 



4i£ {Ei-i£){E2+E3-2i£){E3-3i£) 



.(2) 
"3 

Table 7 

Expressions of the graphs of Fig. 33 The momentum conservation reads A;2h 
A;3+ — for graphs a and /3, /c2+ = + ^3+ for graph 7, and k2+ = ki^ — k. 



ik2 ■ k3)iie){2ks - ki)+ 
ih ■ ks){2ie - E3){ki + ka) 
ih ■ k2)i-E3)i2ki - ks)^ 



3+ 



for graph 6. These expressions must be muhiphed by the overah factor (81). The 



contribution to p_L-broadening is obtained by convoluting with the weight given in 



Eq. (74) 



sign. 

Again, we may put together the graphs with the same kinematics. This is 
shown in Tab. [8] We have left out some imaginary terms, namely 



lm{Ai + Ci-al- /3i) oc 
lm{A2 + C; - as - /Ss) oc 



1 (2fc3 - fci)+ 

ekl(kl{k3+ - ki+) + klh+) ^i+ih - ki)+h+ 

1 (fci + A:3) + 

ekl(kl{k3+ - ki+) + klki+) h+ik^ - ki)+k3+ 



(82) 



which are divergent. The instantaneous-exchange graphs are not enough to 
remove all divergencies, as they do in other cases. But these terms do not play 
any role in our discussion since we must eventually add all complex conjugate 
graphs together. 
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Integration range 


Momentum-dependent factors 












ti 1 


u 1 


Re{Ai + Ci 
Re(^2 + 


— 

— a2 


Si* 

-/32), ^2 


-61 






(2fc3-fci) + 


< A;i+ < ^3+, 
^3+ < h+ < +00 


k\ ■ A;3 


(fci-A:3) + fci+fc3+ 
3(fcl+fc3) + 


^2 C2 p 
1 2 3 


(fci-A:3) + fci+fc3+ 






/33, ^3 - 








(2fci-A:3)+ 




P ^2 p 


(A;i-A;3) + fci+fc3+ 




7! 








k2-kz 


(fcl+2fc3) + 










t2 12 t2 


(fcl+fc3) + fcl+fe3+ 




72 






< A;i+ < +00 


k-i-ki 


3(fcl-fc3) + 








u2 12 u2 

K-^ Kg 

ki-k2 


(fcl+A:3)+fci+fc3+ 
(2fci+fc3) + 




73* 








P2 p u2 


(fc3+fci)+fci+fc3+ 



Table 8 

Sums and differences of the graphs of Fig. 32] and [33] grouped according to kinemat 
ics. The overall constant to be added reads 





ABC 
Fig. 34. px-broadening graphs in which all radiation is in the final state. 
4-2.4 Final state/final state 

We now move on to the case in which all interactions are in the final state 



(Fig. 34). To make the comparison to dipoles easier, we exchange the labeling 
of gluons 1 and 2 in such a way that gluon 2 is always the one whose coupling 
to the quark is closer to the interaction with the target than gluon 1. The 
weights of the graphs A,B,C are shown in Tab. ]9} 

The corresponding dipole graphs are just the complex conjugates of the graphs 



in Fig. 33, i.e. the mirror-symmetric graphs. Thus they do not need to be com- 
puted again and we may directly group the graphs according to the kinematics. 



This is shown in Tab. 10, Note that this time, B has the same kinematics as 7. 



We also see that although the calculation is similar to the case in which all 



interaction vertices are at early times, the resuls shown in Tab. 10 are quite 
more complicated than the ones shown in Tab. |8} 

But again, it is clear that the analytical continuation of S in the ^2+ variable 
is equal to A + C — a — (3 and 5 — 7 in their respective integration domains, 
which proves the equivalence between p^-broadening and dipole amplitudes 
also in this configuration. 
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C2 


< A;2+ < /C3+ 


1 


{ki ■ k3)E2{2k3 - k2) + 


El{E2-E3){Ex+E2-E3) 


C3 






{ki ■ k2)E3{k3 - 2k2) + 



Table 9 

Expressions for the graphs of Fig. 34 where all interactions are in the final state. 



The overall factor is given by Eq. (81) 



Integration range 



Momentum-dependent factors 



Ai + CI 



/3i,i3t-7i* 



(fc2+fc3) 



+ 



^2 + ^2-02-/31, ^2 -72 



< A;2+ < ^3+, 
k-i^ < k2+ < +00 



(k3-k2) + k2+k3 



At + Cs 



On 



1 + 



-1 + 



+ 



(2fc3-A:2)4 



{fc2fc3+)/(fc|fc2+)-l 

4 

l-{kp2+)/{klk3+) 



{kz-k2)+k2+k-i+ 

(ks-2k2) + 
{k-i-k2) + k2+k'j,+ 



< k2+ < +00 



(fc2-fc3) + 
{k3+k2)+k2+k3^ 



{fc2fc3+)/(fc|fc2+) + l 



(2fc3+fc2)4 



-1 + 



l+(A:^fc2+)/(fc^fc3+) 



(fc3+fc2) + fc2+fc3+ 

(fc3+2fc2) + 
(fc3+A:2) + fc2+fc3+ 



Table 10 

Sums and differences of the graphs of Fig.|34|and 33 grouped according to kinemat- 
ics. The transverse momentum-dependent factors are not shown since they are the 
same as in Tab. [s] (third column). Only the real part of the graphs are shown. The 
global constant reads again —QAn'^agNcCp. 



4-3 Four quark-gluon vertices 



We address the case in which there are four quark-gluon vertices in the evo- 
lution of the wave function, linked by two bare gluons. 

Writing the expressions for these graphs is quite straightforward since we treat 

all vertices eikonally: The two gluons both carry the ( ) polarization, and so 

the overall factors need not be discussed since they are identical for all graphs 
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x-x>o<x- 



IIFF (left),IIFF (right) 

Fig. 35. Four quark-gluon vertex graph. This is an interference graph between the 
initial and final state, with 2 couplings in each. There are two possible nontrivial 
cuts. 
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7i 72 
Fig. 36. Dipole graphs corresponding respectively to the left and right cuts of the 



p_L -broadening graph in Fig. 35, As for 72, there are two possible orderings between 



the qg and qg vertices in the complex-conjugate dipole amplitude. 

of similar topologies (except for the minus signs already discussed in several 
occasions which stem from the couplings of the gluons to the antiquark in the 
dipole case). Only the energy denominators are relevant for the comparison. 

We name the p_L-broadening graphs according to the chronology of the four qg 
vertices, distinguishing the graphs in which the loops are nested by a "hat" . 

A relevant way to classify the graphs is to distinguish between interference 
graphs between the initial and final states, and graphs in which the gluons are 
emitted and reabsorbed either in the initial and final state. 



4-3.1 Interference graphs between initial and final state 

Two qg vertices in the initial state, two in the final state. We first 
address the graph in Fig. 35 which has two qg vertices in the initial state, 
and two at late times. There may be either two gluons (leftmost cut) or one 
gluon (rightmost cut) in the final state. These graphs are topologically related 



to the dipole graphs 71 and 72 respectively, shown in Fig. |36| which both 
have 3 dipoles at the time of the interaction. We are going to show that the 
expressions for these graphs are indeed identical. 
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ki+,h -P 




k2+,k2 + k + h 



XX X 



IIFF 

Fig. 37. A p_|_-broadening diagram whose identification with the corresponding 
dipole graph requires analytical continuation. A similar diagram (not drawn) would 
have the cut passing on the rightmost gluon, but the calculation would be the same. 

The energy denominators read 



^iSf E[E2{Ei + E2){E[ + E!2) 
^nF^ ~ ~ E[E'2E2iEi + E2) ~ ~ ^ 



^3) 



The J2 sign recalls that 72 in Fig. [36] actually represents two lightcone pertur- 
bation theory graphs, differing by the ordering of some vertices. The definition 
of the energies may be inferred by applying momentum conservation to the 



graphs in Fig. 35 and 36 The only important point is that the routing of the 
momenta be the same in the p_L-broadening and dipole cases. The "— " sign 



in the second line of Eq. (83) corresponds to the change of one gluon coupling 



from a quark to an antiquark (while for the first graph, two such qg couplings 
become qg couplings, in such a way that the associated minus signs cancel 
each other). 



Note that the first equality in Eq. (83) is straightforward: It is a simple graph- 
to-graph correspondence. We shall not address exhaustively all such trivial 
cases in what follows. 



We now consider the graph shown in Fig. 37 (as well as its dipole partner 



in Fig. 38). We need to write the full expressions of the contribution of the 



graphs to the amplitudes since as we will discover, analytical continuation is 
needed. This is also the opportunity to show how the complete calculation 
goes in the four g (7- vertex case. 

The energy denominators read 



D 



IIFF 



{El - ie){E[ - ie){E[ - E2 - 2ie){E'2 + is) 



B4) 
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where 



El 



2A;i+ 



, E[ 



2h+ 



^2 



2h 



2+ 



, E2 



ih + I2 + hf 



2h 



2+ 



^5) 



(see Fig. 37). The polarization factors for the gluons read 



[d-±ili - p)d±±{-p)d±_{-p - I3)] [d^±{k2)d±±{k2 + h)d^_{k2 + I2 + h)] 



(p-h) ■ (p + h) 



k"^ 



h ■ ih + 12 + h 

k"^ 



The vertices are all eikonal. After having performed the sum over the colors 
of the gluons and averaged over the color of the quark, the associated factor 
reads 



ki+h 



7) 



2+ 



The factors (84), (86) and (87) must be convoluted with the gluon densities as 



72\ 



dki+ dk2+d'^k2 d'^li dH2 d'^l^ xg{x, l\) xg{x, /|) xg{x, I 
37 (27r)3 (27r)3 If If" iV2 - 1 - 1 - 1 



ipLf 



(Momentum conservation was used to get rid of the integration over ki). The 
I factor comes from the ordering of the times at which the gluons interact 
with the nucleus. We may perform the integration over the + component of 
the loop momentum k2 explicitly. The result reads 



dN 
d^p 



Fig. [37] 
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X 



dkij^ 1 
A;f+ {Ei-ie){E[ - is) 



In 



A 



Ck 



1+ 



d^/ca k2 ■ {k2 + I2 + h) 
2vr klik^ + h + hY 



^9) 



The remaining integral over may also be performed. Instead, we stress 
that the content of the square brackets in the last line is the result of the 
integration 



dk2+ . — (90) 

where ( = A;|/(p + ^3)^, and rj = Ae(kl_^_/k2 is a small positive regulator. The 
IR and UV cutoffs are eventually sent to the limits A — )■ and A — >■ +cx). 



We now address the corresponding dipole diagram shown in Fig. 38 We are 
going to show that its evaluation leads to the same expression as in Eq. (89). 
We observe that we may label the transverse momenta of the gluons in such a 
way that only one factor looks different with respect to the previous diagram. 
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Fig. 38. Dipole graph that corresponds to the p_|_-broadening graph of Fig. 37, The 
+ components and /c2+ are always directed to the right. The direction of the 
transverse components is indicated by the arrows above the gluon hnes. There are 
actually two graphs, which are distinguished by the respective ordering of the qg 
and qg vertices in the amplitude (leftmost vertices). Note that this graph is related 
to 72 in Fig. [36] by a simple symmetry. 

namely /+. Indeed, the overall difference between the two graphs is in the 



expression for the energy denominators, which now read (compare to Eq. (84)) 



1 



{El - ie){E[ + ie){E[ + E2 + 2ie)(E^ - ie) 
Then J+ is replaced by 



(91) 



1+ 



I -I- — — / dk2~\- ^ yz , r . 

J\ k2+ {k2+ + Cki+ + ir]) ' 



(92) 



whose explicit expression reads 



/' =ln 



Ck 



— ITT. 



1+ 



(93) 



Hence the dipole amplitude differs from the p^-broadening amplitude by an 
imaginary term, which cancels when one adds the complex conjugate graphs. 

We shall now present an alternative way to view the connection between the 
p_L-broadening diagram and the dipole diagram essentially similar to the an- 
alytical continuation presented before. We remind that this method has the 
advantage that it does not require the explicit evaluation of the integrals, 
hence it is generalizable to cases in which the latter is not possible. It is based 
on the observation that the only difference between the two formulations is a 
change of the direction of the /c2+ momentum. We are going to interpret this 
change as an analytical continuation of the integral /+. We start with the ex- 
pression of /+ in Eq. ( 90 ) . The k2+ plane is represented in the plot of Fig. 39 



The integral in Eq. (90) is represented by the branch along the positive axis 



of the contour in the figure. The Cauchy theorem on the full contour reads 



/+ + /a + I- 



2m 



(94) 
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Fig. 39. Integration contour in the ^2+ plane. The closed contour is divided in 
4 sections: I\ is the small half circle of radius A oriented clockwise, /+ the part 
[A, A] of the real axis, I\ the large half circle of radius A oriented anticlockwise and 
/_ = [—A, — A]. The position of the pole C^i+ + is also shown. 

where the right-hand side is the contribution of the pole at /c2+ = (ki^ + irj. I\ 
is zero when A goes to infinity since the integral is convergent in the ultraviolet, 
and it is easy to see that I\ = in when A goes to zero. Finally, 



dk2+ Ch+ 
A k2+ k2+ - Ch- 



^ dk2+ 

A /C2+ A;2H 



+ Ch+ + iv 



(95) 



Thus, from Eq. (94) and (95), we get 



iTT + lL 



(96) 



which is of course the same equation as (93) and thus leads to the same 
conclusion as before. 



Note that here, we chose to stick to the adiabatic regularization while we 
used the principal part prescription in Sec. 4^ The difference is in imaginary 
terms which anyway cancel when all complex conjugate graphs are taken into 
account. 



Three gluons in the initial state, one in the final state. Let us 

address the case in which there is only one qg vertex after the time of the 
interaction with the nucleus in the p^-broadening case (which corresponds to 
2 real dipoles). The relevant graph (with its two possible cuts) is represented 



in Fig. 40 , and the topologically equivalent dipole graphs are shown in Fig. HT 



We find that the energy denominators for these graphs read 

= -D— = —y-^ = yDs, = y Dg,. (97) 

IIIF IIFI E[El{Ei + E2) ^ ^1 ^ ^2 \ ) 
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IIIF (left cut), IIFI (right cut) 

Fig. 40. p_L -broadening graphs with three qg vertices at early times and one at late 
times, and two gluons (left cut) or one (right cut) in the final state. 
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Fig. 41. Dipole graphs that correspond to the two cuts of the p^-broadening graph 



of Fig. 40 For each of these graphs, the graphs corresponding to different orderings 
of the qg vertices with respect to the qg vertex in the amplitude are understood, 
and this also includes the instantaneous exchange of gluon 2. 

Again, the signs recall that we have included all relevant graphs, see the 
caption of Fig 



41 



The minus sign in front of D — = eventually gets cancelled 
^ iwi_ " 

when one goes to dipoles by the minus sign of the qg vertex in /32. Thus the 
correspondence is verified also for these graphs, namely 



IIIF^^/3i , IIFIm^/32 . 



(9^ 



This is almost a graph-to-graph identity, up to the relative position of some 
vertices in the dipole case. 

We move on to the case in which the loops are disconnected. The relevant 



graphs are shown in Fig. 42 for broadening, and in Fig. 41 for their dipole 
equivalent /Ss and /34. Interestingly enough, there is no graph-by-graph cor- 
respondence although the topologies of IIIF and IIFI seem identical to the 
topologies of /Ss and /34 respectively. 



Indeed, the denominators read 



He[Ex - 3t6){E[ - ie){E2 - le) ' 



(99) 



while D 



ft 



2ie{Ei - 3ie){E[ + ie){E2 - is) ' 
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Fig. 42. Graphs for which 3 quark-gluon vertices are in the initial state, one in the 
final state, and the gluon loops are disjoint. 





/33 



Fig. 43. Dipole graphs that correspond to the p^-broadening graphs of Fig. 42 The 
qg vertex in graph /34 may arrive at any time in ] — oo, 0[; Only the graph in which 
it is posterior to the gluon loop is shown. 

The difference in the global sign is explained by the qg — > qg change when one 
goes to dipoles, but there is an extra sign difference in the energy denominator 
involving E[ which hampers the identification of IIIF and (5^. On the other 
hand, the remaining diagrams have the following denominators: 



D 



IIFI 



2ie{Ei + ie){E[ + ie){E2 - ie) 
1 



2ie{Ei + ie) {E[ - ie) {E2 - ie) ' 
and thus we see that the sum of all graphs satisfy the identity 
_ 1 1 1 

^iiiF + ^iiFi - ,,E^E,^E^ + ^f^^ + = 



(100) 



{Dp, + Y.Dp,). (101) 



Hence the correspondence holds between the sum of the p_L-broadening graphs 
IIIF + IIFl (left-hand side) and the sum of the dipole graphs (^s + J^f^i (right- 
hand side). 



One qg vertex in the initial state, three in the final state. As for 



the case in which there are 3 qg vertices in the final state (Fig. 44), one finds 



similar relations with dipole graphs of the type of those represented in Fig. 45 



and Fig. 46 We may for example check that 



D 



IFFF 



IIIF' 



(102) 
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IFFFJFFFJFFF FIFF,FFIF,FFFI 



IFFFJFFF 



FIFF,FFFI 



Fig. 44. p_L -broadening graphs for the case in which there are 3 qg vertices in the 
final state. 




Fig. 45. Dipole graphs which are topologicahy equivalent to IFFF, IFFF and IFFF 



respectively in Fig. 44 (Similar graphs would correspond to FIFF, FFIF and FFFI). 





Fig. 46. Dipole graphs topologically equivalent to IFFF and IFFF in Fig. 44 
which proves the relation between IFFF and As for the graph IFFF, 

1 



Afff- e[EI{E2-E 



(103) 



which is similar to up to signs, most notably in front of Ei in the de- 

nominator. Analytical continuation in the k2+ variable enables one to identify 
IFFF to IIFL 



D 



IFFF 



i > 



continuation in ^2+ IIFI 



and thus also to using Eq. (97). In the same way, 



D :^ i )■ DxTpT , D ^_ i )■ -Dtttp. 

IFFF continuation in fc2+ FFFI ^^^^ 



The relation with (5^, (5^ then follows from Eq. (99), (100) 



(104) 



(105) 



4.-3.2 No interference between initial and final states 

Finally, we consider the case in which there is no gluon linking an initial-state 
vertex to a final-state vertex. The simplest of these graphs is the leftmost graph 
in Fig. [47l for which a graph-to-graph identification with the leftmost graph 
in Fig. |48| holds (provided one includes all possible time orderings of the qg/qg 
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Fig. 47. Graphs in which the gluons never connect qg vertices in the initial state 
to vertices in the final state. It turns out that the 3 rightmost graphs have to be 
considered together for the correspondence with dipoles to work. 
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Fig. 48. Dipoles graphs that correspond to the p_L-broadening graphs of Fig. 47 As 
usual, this set has to be supplemented by graphs exhausting all possible orderings 
of the vertices, including instantaneous exchanges. The leftmost graph corresponds 
to the leftmost graph in Fig. |47] while the set of the 3 rightmost graphs corresponds 



to the set of the 3 rightmost graphs in Fig. 47 



vertices, which imphes to take into account also instantaneous exchanges in 



Fig. 48 ) . The energy denominators read 

D ' 



nil El{E^ + E2Y 



(106) 



We now turn to a more tricky case, namely the 3 rightmost graphs in Fig. |47 
These 3 graphs correspond to the 3 rightmost graphs in Fig. |48j Indeed, the 
sum of the energy denominators reads 



Ref/^— +£'1111 + AiFF 



1 1 



nn ■ mi. m.y ^2^2 EfE^iE. + E^) (107) 
= -Re^(D«, +D,3 + D„J. 

It may look surprising that there is no graph-by-graph correspondence in this 
case, and especially that one needs to consider together purely initial-state 
graphs such as Illl and a graph where one gluon is in the initial state and 

one in the final state (IIFF). However, we checked that not including the 
latter would lead to a mismatch between p^-broadening and dipoles by a 
term l/{2EfE2). 



We turn to the last class of graphs shown in Fig. 49 This case is solved by 



noting that these graphs may be mapped to the ones shown in Fig. 47, up to 
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Fig. 49. Another set of graphs similar to Fig. 47 in correspondence with the dipole 



graphs in Fig. 48 or with similar graphs. One would also have the graphs in which 
the cut is on the left, symmetric to the rightmost cut. 

an analytical continuation in the variable. We check the identities 
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IFFI continuation in k2-i 
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IIII 



IIII i^FF 



(108) 

and hence the corresponding identities with the (sets of) dipole graphs. Other 
graphs which we do not review here may be deduced from these ones by simple 
symmetries. 

We have left out of the discussion a few graphs which are in exact one-to-one 
correspondence between the two processes. Up to these cases not discussed 
explicitly but trivial, the proof that the evolution with the energy of p±- 
broadening cross sections and of dipole forward amplitudes is identical at 
next-to-leading order is now complete. 



5 Extension to two-jet versus quadrupole amplitudes 



The specific process on which we have focussed so far, namely p^-broadening, 
is the simplest in the class of production processes. However, it is not the 
easiest to measure experimentally. Therefore, we wish to extend our analysis 
to other production processes. 

It has been shown that quadrupole structures appear in the computation of 
observables which involve two particles (dijets) in the final state p^25f26j . 
The energy evolution of quadrupoles can be related to the BFKL evolution 
with saturation, see Ref. |^ (approximation schemes were worked out in 
Ref. [28)129] ). The question is whether the dijet /quadrupole correspondence is 
true when quantum evolution is included. 

The extension of our discussion of the evolution of p_L-broadening/ dipole am- 
plitudes to the evolution of dijet /quadrupole amplitudes is straightforward 
once one notices that the fermion vertices and energies entering the energy 
denominators never play any role. Indeed, the former are represented by a 
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Fig. 50. Two graphs contributing to the evolution of the dijet cross section. The 
incoming object is a qq dipole. 
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Fig. 51. Quadrupole graphs corresponding to the dijet graphs of Fig. [50) On the 
rightmost graph, we have shown how the color lines close to form the quadrupole. 

mere overall factor, which is identical for all graphs of a given class, and the 
latter do not even appear in the calculation, thanks to the soft approximation. 

We illustrate the correspondence between two particular sets of graphs in 



Fig. 50 and 51, which are similar to Fig. 25 and 26 respectively. The relevant 
part of the dijet diagrams may be modeled by an incoming quark-antiquark 



dipole which subsequently evolves by gluon emission (Fig. 50). 



Let us denote by ip the wave function of the dipole, and 5*4 the ^-matrix 
element for the scattering of the quadrupole. Equation ([T]) is to be replaced 
by the following formula (see Fig. SOfsT for the notations): 



dN 



X 5*4(0, X2, a;2, a;'^), (109) 



d?pid?p2 J (27r)2 
which now holds at next-to-leading order in the logarithms of the energy. 



6 Conclusion 



We have shown that the identity ([T]) between the amplitude for p_|_-broadening 
and the dipole cross section holds at next-to-leading order level. The rela- 
tion may be extended to dijet cross sections versus quadrupole amplitudes 
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(Eq. (109)), and probably to other production processes. In other words, we 
have proved that the next-to-leading order BFKL equation describes the en- 
ergy evolution of these production processes. 

Our demonstration was based on a systematic (and quite awkward) inspection 
of all lightcone perturbation theory graphs which contribute to these respective 
processes. We have however avoided the full computation of the graphs since 
our goal was to prove the correspondence in the most general possible way. 

The way the matching occurs turns out to be very subtle. In order to see it, 
we grouped the graphs according to their topologies (number of vertices of 
a given type, chronology of the different interactions) on the p_L-broadening 
side and on the dipole side. The time integrations for the dipole graphs are 
more constrained than for the p^-broadening graphs (there are 4 independent 
integration regions in the latter case, and only 2 in the former case; see the 
discussion of the Keldysh-Schwinger formalism in Sec. [2]). This first apparent 
mismatch is generally solved by subtle cancellations between graphs, and by 
considering appropriate ensembles of graphs (including, in particular, contact 
interactions where gluons are exchanged instantaneously, which turn out to 
play a crucial role). In some cases however, irreducible differences in the flow 
of longitudinal momenta result in differences in the energy denominators that 
can be resolved only by analytical continuation in the relevant longitudinal 
momentum variable. We found a number of pecularities: For example, keeping 
the consistency of the regularization (namely the adiabatic e parameter, see 



Eq. (18)) throughout the calculations turned out to be surprisingly important. 



Despite our efforts, we have not found a simpler argument to explain the 
correspondence more systematically. Such an argument would be crucial for 
example to be able to understand whether the correspondence is still true at 
next-to-next-to-leading order and beyond. 

A full next-to-leading order calculation of production processes requires also 
to take into account the next-to-leading order corrections to the parton dis- 
tribution functions and to the fragmentation functions. In particular, one has 
to release the eikonal approximation at the quark-gluon vertices, which we 
have not done here. Progress in this direction has recently been reported, see 
Ref. IHUfHT] : In these respects, our work may be seen as complementary to the 
latter studies. 
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